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Abstract This paper presents the application of the Szekeres Swiss Cheese 
model to the analysis of observations of the cosmic microwave background 
(CMB) radiation. The impact of inhomogeneous matter distribution on the 
CMB observations is in most cases studied within the linear perturbations 
of the Friedmann model. However, since the density contrast and the Weyl 
curvature within the cosmic structures are large, this issue is worth study- 
ing using another approach. The Szekeres model is an inhomogeneous, non- 
symmetrical and exact solution of the Einstein equations. In this model, 
light propagation and matter evolution can be exactly calculated, without 
such approximations as small amplitude of the density contrast. This allows 
to examine in more realistic manner the contribution of the light propagation 
effect to the measured CMB temperature fluctuations. 

The results of such analysis show that small-scale, non-linear inhomo- 
geneities induce, via Rees-Sciama effect, temperature fluctuations of ampli- 
tude 1CT 7 - 1(T 5 on angular scale i? < 0.24° (t > 750). This is still much 
smaller than the measured temperature fluctuations on this angular scale. 
However, local and uncompensated inhomogeneities can induce temperature 
fluctuations of amplitude as large as 10 -3 , and thus can be responsible the 
low multipoles anomalies observed in the angular CMB power spectrum. 

PACS 98.80.-k, 98.80.Es, 98.65.Dx, 98.65.-r, 04.20.Jb 



1 Introduction 

The Universe, as it is observed, is very inhomogeneous. Among structures 
observed in the Universe are clusters and superclusters of galaxies as well 
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as large cosmic voids. This inhomogeneous matter distribution affects light 
propagation and hence astronomical observations. The study of light prop- 
agation effects is considerably important when analysing the CMB observa- 
tions. This is because the last scattering surface is the most remote region 
which is observable using electromagnetic radiation. In standard approach 
the CMB temperature fluctuations are analysed by solving the Boltzmann 
equation within linear perturbations around the homogeneous and isotropic 
Friedmann-Lemaitre-Robertson-Walker (FLRW) model |T ]l2pl , The applica- 
tion of the FLRW model as a background model results in a remarkably good 
fit to the CMB data [3j. However, the assumption of homogeneity, which is 
also consistent with other types of cosmological observations is not a di- 
rect consequence of them [4]. It is often said that such theorems like the 
Ehlers-Geren-Sachs theorem [5] and the 'almost EGS theorem' [6] justify the 
application of the FLRW models. These theorems state that if anisotropics 
in the cosmic microwave background radiation are small for all fundamen- 
tal observers then locally the Universe is almost spatially homogeneous and 
isotropic. However, as shown in [7], it is possible that the CMB temperature 
fluctuations are small but the Weyl curvature is large. In such a case the ge- 
ometry of the Universe is far from the Robertson- Walker geometry and the 
applicability of FLRW models is not justified. Moreover, the applicability 
of the linear approach can be questionable since the density contrast within 
cosmic structures is much larger than unity. Therefore, there is a need for 
application of exact and inhomogeneous models to the study of the light prop- 
agation and its impact on the CMB temperature fluctuations. This issue has 
been extensively studied within spherically symmetric models — within the 
thin shell approximation [8l l9lfT0] and within the Lemaitre-Tolman model [TT| 
fT ^ fT ^fT4jfT5|fT6ll 1 7| . However, most of the cosmic structures are not spheri- 
cally symmetric, and thus the study of light propagation in non-spherical 
models is essential. One of the suitable models for this purpose is the Szek- 
eres model. The Szekeres model has no symmetry, allows to study a nonlinear 
evolution and does not require small Weyl curvature. Therefore, this paper 
aims to study the CMB temperature fluctuations in the Swiss Cheese Szek- 
eres model. 

The structure of this paper is as follows. Sec. [2] presents basic formulae 
describing light propagation; in Sec. [3] the Szekeres model is presented; Sec. 
|4] presents the construction of models; Sec. [5] presents the findings. 

2 Light propagation 

Light propagates along null geodesies. If k a is a vector tangent to a null 
geodesic, then 

k a k a = 0, k a 0k^ = 0. (1) 

As light propagates, the frequency of photons changes. The ratio of the fre- 
quency of a photon at the emission event to the measured frequency defines 

1 This approac h is im plement ed in such codes like CMBFAST 
(http://www.cfa.harvard.edu/~mzaldarr/CMBFAST/cmbfast.html), CAMB 
(http://www.camb.info/), or (JMBEASY (http://www.cmbeasy.org/). 
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the redshift 

^:=l + z. (2) 

Since photon's energy, as measured by an observer with the 4- velocity u a , is 
proportional to k a u a , thus the redshift obeys the following relation 

where the subscripts e and refer to instants of emission and observation 
respectively. Assuming that the black body spectrum is conserved, the tem- 
perature must be proportional to 1 + z 

£ = ! + (4) 

1 Q 

Then, from eq. (j4]), the temperature fluctuations measured by co-moving 
observer are: 

_ T e /(l + z)-f e /(l+z) 
\t) o f e /(l + z) ' w 

where quantities with bars ~ refer to average quantities. If the temperature 
at the emission is T e = T e + AT e , then 

f AT\ z-z ( AT\ l + z 

As can be seen from the above formula, the observed temperature fluctuations 
on the CMB sky are caused by the light propagation effects and by the 
temperature fluctuations at the decoupling instant. 



3 The Szekeres model 

3.1 The metric of the Szekeres model 

For our purpose it is convenient to use a coordinate system different from 
that in which Szekeres [H| originally found his solution. The metric is of the 
following form [19j 

— /C J hi 

where <P is a function of t and r, e — ±1, and k — k(r) < e is an arbitrary 
function of r. The function E is given by 

E(r,p, q ) = ±(p 2 + q 2 )-^p-Q q + C, (8) 

where the functions S = S(r), P = P(r), Q = Q(r), and C = C(r) satisfy 
the relation 
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Originally, Szekeres considered only a case of p = = A, This result was 
generalised by Szafron [20j to the case of uniform pressure, p — p(t). A spacial 
case of this solution, the cosmological constant, was in detailed discussed by 
Barrow and Stein-Schabes |21j . 

The e = — 1 case is often called the quasihyperbolic Szekeres model (for 
a detailed discussion on the quasihyperbolic Szekeres models see [22]), £ = 
quasiplane (for details see [22153]), and £ = 1 quasispherical (for details see 
[24]). Although it is possible to have within one model quasispherical and 
quasihyperbolic regions separated by the quasiplane regions |22j, only the 
quasispherial case will be considered here. 

In the quasispherial case surfaces of constant t and r are spheres. The 
transformation from (p, q) coordinates into ("&, ip) coordinates is El] 



S cot — cos <p + P, 



S cot — simp + Q. 



(10) 



3.2 The Einstein equations 

Applying metric ([7]) to the Einstein equations, and assuming the energy 
momentum tensor for a dust, the Einstein equations reduce to the following 
two 

?*« = ^-* + ^. (") 

G M, r -3ME, r /E 
* n * P= &{*, r -*E, r /E) - (12) 

where p is matter energy density, M(r) is another function of radial coor- 
dinate. In a Newtonian limit Mc 2 /G is equal to the mass inside the shell 
of radial coordinate r. However, it is not an integrated rest mass but active 
gravitational mass that generates a gravitational field. 
Eq. (fTTjl can be integrated 



^f-k + US 2 



c(t-t B ), (13) 



where ts(r) is an arbitrary function of r. This means that the big bang is 
not a single event as in the FLRW models but occurs at different times for 
different distances from the origin. 

As can be seen the Szekeres model is specified by 6 functions. However, by 
a choice of the appropriate coordinates, the number of independent functions 
can be reduced to 5. 
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3.3 General properties and the Friedmann limit 

The vorticity within the Szekeres model is zero. In addition the acceleration 
vanishes, u";/3 u 13 = 0. The shear tensor is 



1 / $, tr -$, t $, r 



p 3 V ®,r -$E, r /E 
The scalar of expansion is 

fl _ a _ $,tr+2$,t$,r/$-3$, t E, r /E 
" *,r -$E, r IE 



diag(0,2,-l,-l). (14) 



(15) 



The Weyl curvature decomposed into its electric and magnetic part is of 
the following form 

pa r° ,,-y,,' M(3<l>, r —<&M, r /M) (f] 1 n 
£7 p -C lPS u>u - ^—^-^-(11^(0,2,-1,-1), 

- ^^Pcry^C^uV - 0. (16) 

Finally, the 4D and 3D Ricci scalars are 

4 "K = -npc 2 - 4/1, 

, fc ($k, r lk-2$E, r IE \ , 

In the Friedmann limit, <P — ► ra^, fc — > fci?r 2 and M — > Mpr 3 where Of 
is the Friedmann scale factor, /c_f is the curvature index and is a constant, 
and Mp is another constant. As can be seen in Friedmann limit: 

ap ' 



ap 



e - 


+ 3- 






E a (j - 




4 TZ - 




3 TZ - 





3.4 Null geodesic equations 
The geodesic equations 

+ a., 

in the quasispherical Szekeres model, are of the following form 



c 



a = 0: 



dH $, tr -$, t E, r /E 

ds 2 1 — k 

**, t \(dp\ 2 (dq 
E 2 \ ds I V ds 



a = 1: 



(20) 



1 - k 



d 2 r $, tr ~$, t E, r /E dtdr <Z> 

ds 2 <P, r -<PE, r IE ds ds E 2 <P, r -<PE, r /E 



<P, rr -<P, r E, r jE -<f>E rrr jE + $(E, r / E) 2 1 k, 



$, r -$E, r /E 



2 1- k 



<P E, r E, p —EE, pr dr dp £ (E, r E, q -EE, qr ) dr dq 
E 2 $, r -$E, r IE ds ds E 2 $, r -<PE, r /E ds ds 

a = 2: 



(21) 



d 2 p <P, t dt dp 



ds 2 



2—^ C - 

ds ds 



-$E„ IE 



+2 



3V 



£7, 



#(1 - ft) 



dr dp -E,„ / dp 



£ / ds ds 



E 



(E, r E, p —EE, pr ) 



E^dpdq _ 
E ds ds ~ U ' 



E Us. 



a = 3: 



d 2 q t n $, t dt dq 
ds 2 ' 



2-—— — _ 

<P ds ds 



#(1 - fc) 



(£7, r E, q —EE, 



dr dq £7, g / dp 



£ / ds ds E V ds 



2?, p dp dg 



£7 dsds 



(22) 



E, q ( dq 



E \ds 



(23) 



Equations (|20|) - (J23j) are quite complicated. However, if two coordinates 
could be constant on a geodesic, then we could impose dp = dq — on a 
general solution of (7| to get 



dt 
dr 



± 



$, r -$E, r fE 



(24) 



where + is for outwards directed geodesies and — for inwards directed geodesies. 
In such a case the redshift formula (J3j) would reduce to the much simpler form 
(see Appendix [A] for derivation) 



ln(l + z) = ± d 



$, ir -$, t E, r IE 

VT^k 



(25) 



Thus, if such fixed-direction geodesies exist, then the study of light propaga- 
tion in the Szekeres model could be significantly simplified. Instead of solving 
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eqs. |20ll - (|23|) only eq. |24|) would have to be solved to find a null geodesic 
and only eq. ([25| to find the redshift. However, in general d 2 p/ds 2 ^ and 
d 2 gr/ds 2 7^ 0, i.e. the condition dp = dq = cannot hold along the null 
geodesic. As follows from ([20)1 - |23j) if initially k p — k q — 0, then the coor- 
dinates p and q will remain constant only if along the whole geodesic 

<*V = <^, (26) 

or 

EE, pr — E, p E, r and EE, qr = E, q E, r . (27) 

The relation lj26|) holds only at a shell crossing singularity which must be 
eliminated in a physically acceptable model. Apart from the spherical sym- 
metry (i.e. the Lemaitre-Tolman model) the relations ([2"7[) hold only in the 
axially symmetric case [25j. It should be noted that not every Szekeres model 
is axially symmetric. Moreover, apart from the spherical symmetry there is 
only one such geodesic, which propagates along the symmetry axis. Thus, 
such a geodesic will be referred to as the axial geodesies. 



4 The Swiss Cheese model 

4.1 Arrangement of the Swiss Cheese model 

The Swiss Cheese models which are employed in this paper are constructed 
from six different building blocks - regions A-F (holes) - which are matched 
with the Friedmann background (cheese). These Szekeres patches are placed 
so that their boundaries touch wherever a light ray exits one inhomogeneous 
patch. Thus the ray immediately enters another Szekeres inhomogeneity and 
does not spend any time in the Friedmann background. Using different se- 
quences of regions A-F five models are constructed, . The density distribu- 
tion at the current instant within each of these regions is presented in Fig. 
[TJ The exact forms of functions used to specify the Szekeres model in each of 
these regions are presented in Appendix [B] As can be seen from the detailed 
specification (Appendix [Bj the functions defining regions A-D become for 
the radial cooridat^l r ^ 24 kpc of exactly the same form as the form in 
the Friedmann background [compare the form of the functions in Appendix 
[B]with the form of functions in the Friedmann limit, eqs. lfl8|) ]. Regions E 
and F tend exponentially to the Friedmann models. However, as seen from 
their specification or from Fig. [2j at the distance r ~ 30 kpc and r w 40 kpc, 
respectively, regions E and F become almost Friedmann. Figure [2] presents 
the curvature scalar, W 2 , which is defined as 

W 2 -^g^, (28) 

2 The radial coordinate in this paper is defined by the value of <P at the last 
scattering instant: r = <1>ls - see Appendix[B] Thus r = 24 kpc corresponds to the 
current distance of approximetaly 26 Mpc - cf. Fig. [l] 
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Fig. 1 The density distribution, p/pb, at the current instant for regions A (up- 
per left), B (upper right), C (middle left), D (middle right), E (lower left) and F 
(lower right). In regions A, B, E and F the dipole axis is aligned with the Y axis. 
For these regions the density distribution is presented for the surface Z = 0. In 
regions C and D the dipole is aligned with the Z axis. For these regions the density 
distribution is presented for the surface Y = 0. Coordinates X, Y, Z are denned as 
X = ^>cos Lp sin 1?, Y = <£sin </Psin Z = ^cosi? [where $ and ip are given by l|10p ]. 
The exact forms of functions used to specify these models are included in Appendix 

m 
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where E a p is the electric part of the Weyl tensor (fl6l) and H = (1/3)8 is 
the Hubble parameter lfl5|) . As can be seen, in some regions W 2 ^> 1. This 
feature, apart from nonlinear evolution and non-symmetrical shape, makes 
the application of the Szekeres model more realistic. 



4.2 Junction conditions 

When constructing a Swiss Cheese model, we need to satisfy the junction 
conditions for matching the particular inhomogeneous patches to the Fried- 
mann background, and also assure the continuity of the null geodesies. The 
standard junction conditions are that the 3D-metric of the surface and its 
extrinsic curvature, the first and second fundamental forms, must be contin- 
uous. 

For matching a Szekeres patch to a Friedmann background across a co- 
moving spherical surface, r — constant, the conditions are: that the mass 
inside the junction surface in the Szekeres patch is equal to the mass that 
would be inside that surface in the homogeneous background; that the spatial 
curvature at the junction surface is the same in both the Szekeres and Fried- 
mann models, which implies that k$z — kpr 2 and (ksz),r = Skpr; finally 
that the bang time and also A must be continuous across the junction. The 
mass M and the curvature function k are matched by the construction — see 
Appendix[Bj The value of the cosmological constant is the same in all regions, 
and the value of the bang time function is fixed by (fl3l) . and at the junction 
is equal to ts = 0. It might be surprising that a non-symmetrical model like 
the Szekeres can be joined with the symmetric FLRW model, but there are 
other examples of such junctions. For example Bonnor demonstrated that 
the Szekeres model can be matched to the Schwarzschild solution [26J. 

The junction of null geodesies requires the continuity of all components 
of the null vector. However, let us notice that when one Szekeres sphere is 
matched to another Szekeres sphere it can be rotated around the normal 
direction. Thus, we only need to match up the time component, k° and the 
tangential component. The tangential component is defined as 

k T = %^(kP) 2 + (ki) 2 . (29) 
hi 

The radial component is then given by the null condition, k a k a = 0. 



4.3 Description of models 

Five different Szekeres Swiss Cheeses models are considered here: 
1. Model 1 

Model 1 is constructed from alternately matching regions A and B (A + 
B + A + B ...) into the Friedmann background. When a light ray exits 
one Szekeres region, it immediately enters another inhomogeneous patch. 
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Fig. 2 The Weyl curvature scalar W 2 [eq. (j 1 6 j> ] evaluated at the current instant, 
along the dipole axis, R, in Szekeres regions A-F (for regions A, B, D, and E, R 
=Y, for regions C and D, R=X). The 3D shape of these profiles is very similar 
to the shape of the density distribution - see Fig. [T] Since the FLRW models are 
conformally flat, the Szekeres regions are in some parts far from being even a close 
(linear) approximation of the FLRW model. 



Each time the (p, q) position of the point of entry is randomly selected. 
In addition k p and k q are quasi-randomly selected, i.e 

( fc9 ) 2 =^ T f) 2 < ( fcP ) 2 = (l-7)(fc T f) 2 , 

where 7 is a random value in the range 0^7^ 1. The radial coordinate 
of the matching point is 77 — 24 kpc - the point where the Szekeres region 
becomes Friedmann. 

2. Model 2 

This model is constructed from alternating regions C and D, but only 
axial null geodesies are considered, i.e. k p = and k q = 0, p = q = 0. 
The radial component of the matching point is again 77 = 24 kpc. 

3. Model 3 

The next model consists of regions E and F placed alternately. Null vector 
components k p and k q are chosen in such a way that 10~ 8 ^ k p ^ 10~ 4 
and 10~ 8 < k q ^ 10~ 4 , but are otherwise random. As can be noted, this is 
not in accordance with condition lj29j) . In order to maintain the continuity 
of the tangential component of the null vector the next Szekeres patch 
must by reoriented with respect to the preceding patch. This however 
leads to an overlapping of these two Szekeres regions. Although, at the 
junction point (77 — 40 kpc for region E, and 77 = 50 kpc for region 
F), these two regions are almost Friedmann, still this is not a perfect 
matching. We proceed with this type of imperfect matching to study 
how randomly chosen values of the tangential component (hence more 
randomised light propagation through a Szekeres patch) influences the 
final results. 

4. Model 4 

Model 4 is constructed using only C regions, with 77 = 24 kpc, and only 
axial geodesies are considered, i.e. k p — k q = p = q = 0. 
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Fig. 3 The temperature fluctuations caused by light propagation effects in models 
1-5. In models 1-3 light propagates alternatively though underdense and overdense 
regions. In model 4 light propagates only though regions of 5M < 0, 5k > 0, and in 
model 5 only through regions of 5M > 0, 8k < (see Sec.|4]for a detail description 
on how these models were constructed). 




Fig. 4 A small part of the light propagation in model 3. The left panel shows 
the density variation that light 'feels' as it propagates. The black thin dotted line 
shows the density in the background model. The right panel presents the tem- 
perature fluctuations as measured by an observer situated outside the structure 
in the homogeneous FLRW region. The letters in the left and right panels label 
corresponding points along the light path. 



5. Model 5 

The last model is also axially symmetric, k p = k q = p = q = 0, rj = 24 
kpc, but uses only D regions. 



12 



5 Results 

5.1 The Rees-Sciama effect 

To estimate the temperature fluctuations induced by the light propaga- 
tion effects, it is assumed that initial temperature distribution is uniform, 
(AT/T) e = 0. Then temperature fluctuations are calculated using (??), and 
they are plotted against time of propagation in Fig. 03 As seen, the final val- 
ues are small, of amplitude AT jT « 1CT 7 (model 3), AT jT « 1(T 6 (models 
1 and 2), and AT jT « 10~ 5 (models 4 and 5). A detailed analysis of how 
inhomogeneities induce temperature fluctuations is presented in Fig. [4] (for 
clarity, only a small fraction of the time is presented). The left panel of Fig. 
|4] shows the density of regions through which the light propagates in model 
3. The right panel presents the temperature fluctuations as measured by an 
observer situated at the junction point where the model is that of Friedmann. 
Letters correspond to each inhomogeneous patch (left panel) and tempera- 
ture fluctuations caused by them (right panel). Clearly, underdense regions 
induce negative temperature fluctuations, overdense regions induce positive 
fluctuations. 

Apart from estimating the amplitude of the Rees-Sciama effect, it is also 
important to estimate the angular scale which is the most affected by this 
effect. Without going into any complicated analysis, we can estimate the 
angular scale by employing the following approximation: the correlation be- 
tween two distant points on the sky is zero - photons which were propagating 
along two distant paths have the temperature fluctuations uncorrelated. Only 
when the light paths are near to each other are the temperature fluctuations 
correlated. Thus the simplest estimation of the angular scale of the Rees- 
Sciama effect, as seen from the schematic Fig. [5j is the angular size of the 
Szekeres patch at the last scattering instant. For the models studied in this 
section, such approximations lead to an angular scale of i) 0.21°, or alter- 
natively I « 850. If the photons are propagating along neighbouring paths 
for only half of the age of the Universe (in such a case, as seen from Fig. 
[BJ the final temperature fluctuations are two times smaller), then the angu- 
lar scale is similar, d s» 0.24° (£ ss 750). Thus, the Rees-Sciama effect of 
amplitude ~ 10~ 6 contributes to the CMB temperature fluctuations on the 
angular scale ?9 < 0.25° (£ > 700). This angular scale corresponds to the 
angular scale at which the third peak of the CMB angular power spectrum 
is observed. At this scale the measured rms temperature fluctuations are of 
amplitude «5x 10 -5 . This is still several times higher than the results ob- 
tained within models 4 and 5. In the case of models 1-3 the measured value 
is more than one order of magnitude larger than the model estimates. 

5.2 The role of local structures 

So far, it has been assumed that each inhomogeneous structure is compen- 
sated (i.e. each Szekeres region was matched with the FLRW background), 
and that measurements are carried out away from the inhomogeneities, i.e. 
where the universe is homogeneous. However in the real Universe there is no 
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last scattering surface 




now 

Fig. 5 The schematic representation of the Swiss Cheese model. When two photons 
are propagating along a similar path the final temperature fluctuations are similar. 
If paths are different, then the final temperature fluctuations are also different and 
hence not correlated. 
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Fig. 6 Temperature fluctuation amplitude, as measured by observers at different 
locations in region E, along the path of a light ray. The ray path is shown as the 
bright line in the upper inset. 
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place where the cosmic structures in the observer's vicinity are fully com- 
pensated and therefore the Universe should not be treated locally as ho- 
mogeneous. Since all measurements are always local, let us consider what 
happens if temperature fluctuations are measured in an uncompensated re- 
gion. Figure[6]presents the temperature fluctuations measured by an observer 
situated at different places within region E. These results are obtained under 
the assumption that light from last scattering is propagating through ho- 
mogeneous regions, and currently reaches an observer in an inhomogeneous 
structure (region E). The light enters and propagates along the bright line 
shown in the upper right inset of Fig. El The above results show that local 
structures can significantly contribute to the CMB temperature fluctuations. 
This indicates that care must be exercised when extracting information from 
the CMB observations. Although it is highly unlikely that the signal caused 
by the local structures have a signature of acoustic oscillations we should be 
aware that local structures can have some visible impact on observations. 
Thus, it is important to test if local structures can cause the observed corre- 
lations of the alignment of dipole, quadrupole and octopole axes of angular 
power spectrum of the CMB temperature fluctuations (cf. [27.28,29,3011311 
[32] ) or their low amplitude (cf. [33" ll34P35| ). 



6 Conclusions 

The analysis presented in this paper has aimed to examine the influence of 
the light propagation effects on the temperature fluctuations of the CMB. 
The results indicate that the Rees-Sciama effect caused by the propagation 
of light through inhomogeneous but compensated structures do not signifi- 
cantly affect the CMB temperature fluctuations. Some would say that this is 
an obvious result since similar conclusion is reached when using perturbative 
methods. However, as it was argued in the Sec. 03 within real cosmic struc- 
tures the density contrast and the Weyl curvature are significantly large. 
Thus, the application of the perturbation methods cannot be justified. It 
was shown, in this paper, that even in such cases light propagation effects 
are smalH. It has also been shown that the Rees-Sciama effect of amplitude 
~ 10~ 6 contributes to the CMB temperature fluctuations on the angular 
scale i9 < 0.25° (£ > 700). 

However, if the structures are not compensated or the measurements are 
carried out inside inhomogeneous non-compensated structures, the amplitude 
of measured temperature fluctuations can be slightly higher. Since in reality 
we cannot separate ourselves from the surroundings and say that all local 
structures at our positing are "compensated", thus the local cosmic structures 
must be taken into account when analysing CMB observations. Especially, it 
is possible that the local structures can have some impact on low multipoles 
anomalies of the angular CMB power spectrum. 



3 However, in this paper all inhomogeneities are at the current instant of diame- 
ters of ~ 50 — 60 Mpc. Larger inhom ogeneities, of diameters ~ 600 Mpc, can have 
more sig nificant impact, cf. Pl0| ll7|. 
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A The redshift formula for axial geodesies 

To find a simplified redshift formula for the axial geodesies first we have to find k° 
in an affine parametrisation and then use the relation J3]). Since we assume that 
k 2 = fc 3 = then let us choose 

VT^k y ' 

The above parametrisation is not affine, thus the parallel transport does not pre- 
serve the tangent vector. Therefore k°", after being parallely transported, becomes 
Xk a (where A is a scalar coefficient - a function of the parameter s along the 
geodesic). In this case the geodesic equations are of form [36J 

k a ; k? = -~k a , (31) 

and reduce to 

1 dA $, rr -$, r E, r IE - $E,„ IE + ${E, r / E) 2 



A As &, r -$E, r IE 

, n $,tr-<P,t E, r /E k, r 

All the quantities above are evaluated on the null geodesic. In order to better depict 
which quantity is evaluated on the null geodesic a symbol * is be used. Since on 
the geodesic t and r are connected with each other via relation lf24|) . we have 

$n,r = )n"p + (#,r )n 

{$,r ) n ,r = ($,tr )n4~ + (*,rr )n, (33) 

where the subscript n referees to qua ntities measured on the geodesic. 

The second term in equation ll32|l looks like a logarithmic derivative. However 
because of (1 3 3 f) we have 

din [{$, r )n-$n{E, r / E) n ] = 

dr 

(0, r ) n , r -0 n , r (E , r / E) n - $ n {E, rr / E) n + $ n (E , r / E)j = 
(<P, r ) n -<P n (E, r /E) n 

($,rr )n ~ (#,r )n(E, r / E) n - $ n (E, rr / E) n + $ n (E, r / E) 2 n 



($, r ) n -$ n (E, T /E) n 
(<P,tr)n ~ (g,t)n(£,r / E) n At 
{p,r)n-$n(E, T /E) n Ar' 



(34) 
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Using the above relation we can integrate equation ([32]) 



Now we can easily find that k a in the affine parametrisation is given by k a = 
{\/C)k a . Using © we obtain 

ln(l + z)=± T L r ^) n -{^) n {E, T /E) n 
J VI - k n 

r s 

where + is for r e < r and — for r e > r . Alternatively the redshift can be found 
by integration over time: 



B Model specification and evaluation 

In order to define the Szekeres model five functions of radial coordinate needed 
to be specified. In this paper all models will be defined by the following set of 
functions: k, M, S, P, and Q. 

The algorithm used in the calculations can be defined as follows: 

1. The radial coordinate is chosen to be the areal radius at the last scattering 
instant r = <£(r, tLs)- However, for clarity in further use, the prim is omitted 
and the new radial coordinate will be referred to as r. 

2. The chosen background model is the ACDM model, i.e. a flat FLRW model 
with A 7^ 0. The background density at the current instant is then given by 

p b = Q m x p cr = 0.27 x (38) 

07TLr 

where the Hubble constant is f/o = 72 km s^ 1 Mpc" 1 . The cosmological con- 
stant, A, corresponds to Qa = 0.73, where Qa = (l/3)(c 2 A/Hq). 

3. The initial time, ti, is chosen to be the time of last scattering, and is calculated 
from the following formula for a background FLRW universe [37] 

oo 

t(z) = —j- f dz - , (39) 

Ho J (l + z)^V(S) 

z 

where: 

V(z) = Q m {\ + zf + Q K {1 + zf + Q A , (40) 

where flu = 1 — Om — ^A- For the lower limit of integration, z = 1089 was 
used as the redshift at last scattering. 

4. Six different Szekeres regions are considered in this paper. Let us denote them 
as regions A, B, C, D, E and F. The functions M, k, Q, P and S in these 
regions are defined as follows 

— regions A and B 

Mir 3 for r sC 0.5a, 

w = u, - i M * ex P [- 12 (^) 2 ] for °' 5a < r < L5a 
Mi (2a - r) 3 for 1.5a ^ r ^ 2a, 

, for r ^ 2a, 
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where M b is the mass in the corresponding volume of the homogeneous uni- 
verse, i.e. M b = (47rG/3c 2 )p LS r 3 , p LS = pb{l + z LS )'\ Mi = 8M 2 a" 3 e- 3/2 , 
M2 is equal to —0.3 kpc and 0.2 kpc for region A and B respectively, and 
a = 12 kpc. 



( fcir 2 for r < 0.56, 



1 ; 

k = x < 

2 



k 2 exp [-4 (^) 2 ] for 0.56 sC r < 1.56 

fci(26-r) 2 for 1.56 < r < 26, 

I for r ^ 26, 

where fci = 4fc 2 a -2 e -1 , fc 2 is equal to —5.15 x 10 -6 and 3.5 x 10 -6 for 
regions A and B respectively, and 6 = 10.9 kpc. 

5=1, P = 0, Q = Q 1 ln(l + Q 2 r)xexp(-Q 3 r), 

where, for regions A and B respectively, Q\ equals —0.72 and —1.45, Q 2 
equals 1 kpc -1 and 0.4 kpc -1 , and Q3 equals 0.01 kpc -1 and 0.005 kpc -1 . 
With these definitions the mass distribution and the curvature are the same 
as in Priedmann models, for r > 24 kpc. 
Regions Ci and C2 

In region C the functions M and k are the same as in region A. The only 
difference is in the form of functions 5, P, and Q which are as follows 

S = e ar , P = Q, Q = 0, 

where a equal to —0.0255 kpc -1 and +0.0255 kpc -1 for regions Ci and C 2 
respectively. Region Ci is the mirror image of C2, where the Z = surface 
is the symmetry plane [Z = <?cosi9 and $ is defined by the stereographic 
projection fllOp j. The reason for employing two mirror-similar regions is 
that in the coordinates used here, the axial geodesies can only be studied 
for propagation along the Z < direction, in which $ = — ir. Along the 
Z > direction we have $ = 0, which corresponds to a point at infinity 
in the stereographic projection. This problem is overcome by matching Ci 
with C2 along the surface of Z = 0. When calculating propagation toward 
the origin model Ci is employed, and when calculating propagation away 
from the origin model C2 is employed. In both models light propagates 
along the Z < axis. 
Regions Di and D 2 

In region D the functions M and k are the same as in region B. The only 
difference is in the form of the functions S, P, and Q which are of the 
following form: 

S = r a , P = 0, Q = 0, 

where a equal to —0.97 and +0.97 for regions Di and D2 respectively. As 
above, region D comes from matching regions Di and D2 along the Z — 
surface. 



Regions E and F 

Mir 3 for r < 0.5a, 

M 2 exp [-6 (^) 2 ] for r > 0.5a, 



M = M b + 



k = 



kir 2 for r ^ 0.56, 



exp [-(£+!) 2 ] for 0.56, 

5=1, P = 0, Q = Qi-0.221n(l + Q 2 r)xexp(-Q 3 r). 

where Mi = 8a -3 M 2 e -1 ' 5 , ki = 4a -2 fc 2 e -1 . For region E, M 2 = -0.75 
kpc, a = 15.23 kpc, k 2 = -1.00173 x 10 -5 , 6 = 12.95 kpc, Qi = -0.22, 
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Q 2 = 1 kpc -1 , Q 3 = 0.1 kpc -1 . For region F, M 2 = 0.9 kpc, a = 23.76 
kpc, k 2 = 7 x 10~ 6 , b = 19.1 kpc, Qi = -1.4, Q 2 = 0.4 kpc -1 , Q 3 = 0.005 
kpc -1 . 

5. Light propagation was calculated by solving eqs. ll20l - ll23l (models 1 and 3) 
and l|24p (models 2, 4 and 5) simultaneously with the evolution equation l|lip . 
At each step the null condition, k a k a = was used to to test the precision 
of calculations. All equations were solved using the fourth order Runge-Kutta 
method. 

6. The temperature fluctuations were calculated from eg. (J6j) . The redshift was 
calculated using relation ((3J (models 1 and 3) and l]2"i[| (models 2, 4 and 5). 
The mean redshift z was calculated using the ACDM model. 
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